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A closed expression for the momentum evolution of a test particle in weakly-coupled plasma is 
derived, starting from quantum many particle theory. The particle scatters from charge fluctuations 
in the plasma rather than in a sequence of independent binary collisions. Contrary to general belief, 
Bohr's (rather than Bethe's) Coulomb logarithm is the relevant one in most plasma applications. A 
power-law tail in the distribution function is confirmed by molecular dynamics simulation. 
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Though Coulomb scattering is a most basic process in 
plasma and has been studied for a century yLj, doubts 
concerning the treatment as a sequence of independent 
binary collisions remained 0, 13 , and recent analysis 
has revealed that this standard assumption is not 
justified in general and requires revision. Here we de- 
rive the time-dependent many-particle wavefunction of 
a test particle simultaneously interacting with N parti- 
cles residing in the Debye sphere. The plasma parame- 
ter N = nX 3 > 1 involves the plasma density n and the 
screening length A = max(vT /w p , Vq/u> p ), where vt is the 
thermal velocity of ele ctrons, vq the velocity of the test 
particle, and lo p — \J Airne 2 / m e the plasma frequency. 
We emphasize that the collective interaction described 
here is different and in addition to Debye screening; it is 
not included in the usual dielectric approach [j| . The new 
results can be viewed as interaction of the test particle 
with the charge fluctuations inside the Debye sphere; in 
this picture, the test particle of charge Zoe is scattered by 
an effective, spatially extended charge ey/~N rather than 
by a sequence of binary collisions. 

This has deep consequences for the Coulomb loga- 
rithm, because it drastically shifts the borderline be- 
tween classical and quantum Coulomb scattering, ex- 
tending the domain in which the classical approxima- 
tion applies. This is shown in Fig. ^ In the binary 
collision approach (Fig. [IJi), the borderline is given by 
the parameter a = Zae 2 /hvo such that a < 1 defines 
the quantum-mechanical region where the Born approx- 
imation leads to Bethe's logarithm L q = ln(Xm e Vo/}i) 
0, while for a > 1 classical mechanics apply leading to 
Bohr's logarithm L c i — ln(Am e «Q / Z$e 2 ) . In the present 
theory instead, the borderline is given by aN 1 / 2 f» 1, 
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FIG. 1: Regions in a density-temperature plane (atomic 
units) in which Bohr's classical Coulomb logarithm (white 
area) and Bethe's quantum expression (hatched area) ap- 
ply; (a) binary collision theory with borderline defined by 
a — Zoe 2 /hvo = 1, and (b) present theory with approximate 
separation along a 2 N{l + In AT) = 1. Also shown as grey 
area are the region of strongly non-ideal plasma (borderline: 
T ~ n 1/3 ) and the region of degenerate plasma (borderline: 
T~n 2 / 3 ). 



and this leads to a very different picture in Fig. lb. Now 
L c i applies to almost the entire high-temperature region, 
including e.g. the important domain of magnetic fusion 
plasmas, while L q plays only a marginal role. 

Let us first discuss this result in qualitative terms. A 
particular feature of the Coulomb (1/r) interaction is 
of crucial importance in this case, namely that scatter- 
ing does not depend on a and ft, as we know from the 
Rutherford cross section. The difference between a < 1 
and a > 1 regions arises only when the potential devi- 
ates from 1/r, as it is the case in a plasma due to De- 
bye screening occurring at long distances A. This means 
the distinction between classical and quantum treatment 
reveals itself for small-angle scattering, while close colli- 
sions with large-angle scatter are not affected. This point 
has been emphasized strongly by Bohr (see p. 448 in [l|) 
and also in Sivukhin's review [3( (p. 109-113). In accor- 
dance with Bohr "any attempt to attribute the difference 
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between [the classical a ^> 1 and quantum a <C 1 cases] 
to the obvious failure of [the classical] pictures in ac- 
counting for collisions with an impact parameter smaller 
than [the de Broglie wave-length] will be entirely irrel- 
evant. In fact, this argument would imply a difference 
between two distribution for the large angle scattering, 
while the actual differences occur only in the limits of 
small angles." 

Now let us compare the classical scattering angle 5 c i = 
(Z ee*/moV 2 )/\ at distance A with that of quantum 
diffraction S q = (h/m v )/X 0. The open question here 
concerns the effective net charge e* which the test particle 
experiences when passing the Debye sphere. The value of 
e* is not evident, because we are dealing with Coulomb 
collisions at distances much larger than 1/n 1 / 3 . The bi- 
nary collision approximation circumvents this predica- 
ment by alleging that the total scattering can be treated 
as the sum of independent binary interactions happening 
at different times One is then led to take e* = e for 
granted instead of actually calculating e*. The central 
result of this paper will be that the effective charge is 
essentially given by « eN 1 / 2 . The matching condition 
then is S c i/S q « aN 1 / 2 w 1 replacing the condition awl 
obtained in binary collision approximation [3|. The the- 
ory underlying Fig. lb will now be derived. As a central 
result, we also present molecular dynamics simulations 
confirming the analytic theory. 

The present analysis starts from a full quantum- 
mechanical description of the plasma in terms of the 
many-particle wave-function %p — exp(iS/h). The action 
function S satisfies the equation 

dS ^({VjS) 2 ^ TT AjS\ 

— ^=> + > u ik -ih—?—\, (l) 

dt ^ I 2m, ^ °' k 2m, I ' V ; 

3 \ 3 k>j J J 

where the indices j and k denote plasma particles for 
j, k = 1,2,3... and the test particle for j, k = 0, rrij are 
the masses, and Uj tk represent the Coulomb interactions. 
Eq. is equivalent to the exact Schrodinger equation 
and has the form of a Hamilton- Jacoby equation with 
additional terms that are proportional to H and describe 
quantum effects. We examine the solution of Eq. Q for 
the particular initial conditions S(t — 0) = X]j>o Pj ' r i' 
(the Green function for the coordinate-momentum rep- 
resentation), where Pj are given momenta of the plasma 
particles at t = 0, and introduce 

a = S - (Pi • Tj - p)tj2m 3 ) 

such that a(t = 0) = 0. Now the clue for solving Eq.QJ 
is the high-energy approximation 7] which applies to an 
almost ideal plasma and requires |pj| 3> |Vj<t|. Under 
this approximation we find 

j \ k>j J I 



where Vj = pj/rrij. The solution of Eq. (01 for a(t — 
0) = is 

CT = -EE f Ukj (T>kj) dr, (3) 

j>0k>j J ° 

where Sr ktj = r k - v v 5v ktj = v k - Vj and B kJ = 
Sr k j — 5~v k j (t — r). It can be verified by direct insertion. 
The terms Aja proportional to h vanish for the special 
case of Coulomb U oc 1/|D| interactions for distances 
|D|>0. 

The problem of solution @ is that it contains the sin- 
gularities of U k ,j (Dfcj ) for close encounters with \D k ,j \ — > 
0. This deficiency is due to the high energy approxima- 
tion. Inserting Eq. © into \pj\ 3> |Vj<t|, we find that 
this condition is fulfilled only for regions 

\T>k,j\ > \eje k \/(p kj \Sv kj \ 2 ) (4) 

with < t < t and ji k ,j — mjm k /(mj + m k ). Had 
we solved the nonlinear equation Q exactly, we had 
obtained a non-singular result with the maximum mo- 
mentum transfer of 25vj tk fJ'j,k, a s we know from Ruther- 
ford scattering. The way to deal with this problem 
is to cut out in the wavefunction those spatial regions 
which do not satisfy Eq. J3J. The cut-off @ warrants 
that the maximum momentum transfer 2Svj^ k fj,j tk is pre- 
served; this can be verified by operating with —ih\7j 
on exp(iS/H). It should be understood that this short- 
range cut-off is a technical correction (compare with 0, 
p. 448-449:". . . the central region of the field. . . , which, 
on classical mechanics, is responsible for all large angle 
scattering will, for a <C 1,... gives rise to only a frac- 
tion of the order a 4 of the Rutherford scattering"). It 
has nothing in common with differences between clas- 
sical and quantum scattering. These reveal themselves 
only at long ranges 0, ■ Another detail concerning the 
general wavefunction concerns the initial conditions. In 
case the initial state of the plasma is defined by the wave- 
function <f)(t — 0, pi, P2, . . • ) rather than by a fixed set 
of momenta, the corresponding general wavefunction is 
given by 

-0(i,r o ,ri, . . . ) = / exp(iSy/i)(Ki = 0,pi,p2,...)IJ dp k . 

J k>l 

As it turns out, the explicit form of tfi is of no relevance 
in the applications discussed below. 

We now have at our disposal in analytical form the 
time-dependent many-particle wavefunction describing a 
dilute high-temperature plasma. This is a remarkable 
result. An outstanding feature is that it describes simul- 
taneous multiple Coulomb interactions between the par- 
ticles and, in this respect, goes beyond the binary colli- 
sion approximation. Another essential property is that it 
holds for both the quasi-classical regime (a ^> K) as well 
as the deeply quantum-mechanical regime (a <ti h) and 
therefore provides a unique tool to investigate the transi- 
tional region. We now proceed to use this wavefunction 
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to calculate plasma properties. This is straightforward, 
though tedious, and therefore we can give here only the 
main results, leaving technical derivations to a separate 
publication. 

We first consider the distribution function M(t, Q) of 
transverse momentum Q of a test particle moving at time 
t = with momentum p collinear to the x-axis. For 
brevity we consider a fast ion ttiq 3> m e , vq 3> v t for 
times that are longer than to — X/vq, though shorter 
than the collision time t c ~ moNto/rn e L c i, q . M(t, Q) is 
obtained as the matrix element 



M{t, Q) = J exp (iQ • R/K) F(t, R) d 2 R/(2Trhf 



(5) 



where 



F(t, R) = ~ f ^(t,r ,r 1 ,...)^*(«,r +R,ri,...)n ^ 

J fc>0 

and V is the plasma volume. 

Expression |J5} can be significantly simplified for the 
case under consideration. The test particle affects only 
plasma particles in the interaction sphere \Stqj \ < A, for £ Qr q2 < anc j 
which two-particle correlations among plasma particles 
are small owing to T > e 2 ^ 1 / 3 . Aiming deliberately for 
calculations with logarithmic accuracy, we can omit the 
integration over |£roj-| > A and use the method devel- 
oped in 0. We find F(t, R) = F e (t, H)Fi(t, R) where 
F e (t,TL) = exp(-/ e (i,R)), 



quantum regime is restricted to 7 < 1, while the clas- 
sical regime is found for 7 > 1 and will be discussed 
first. Evaluating Eq. © in the limit of very small R, 
one finds / e (R) = vt ln(Am e w 2 /Zoe 2 )(R//i) 2 with v = 
2irnZ 2 e 4: /vo. This is the relevant region in the Fourier 
integral of M e for large enough time t (t\ -C t < t c , 
see below). The function F e (R) is then a Gaussian, and 
M e can be easily calculated. Setting Fi = 1, we obtain 
(Q 2 ) e = J Q 2 M e {t, Q) d 2 Q = 4vtL cl and recover the 
classical Coulomb logarithm L c i = ln(Xm e VQ / Zge 2 ) . The 
important new result here is that it applies to the whole 
region 7 > 1 and not just to a > 1. 

It should be noticed, however, that the function 
/ e (R) = vt ln(A/i?)(R/fi) 2 is more complicated in gen- 
eral and contains a factor lni? for larger radii, where 
R = max(a\R\,ah/m e vo) for a > 1 and R = 
max(|R|, ah/m e vo) for a < 1. For short times, just some- 
what larger than X/vq, this logarithmic factor modifies 
the Gaussian character of M e (t, Q), giving it a power- 
law tail at high Q = |Q|. We then obtain 



M(*,Q)=exp(-Q 2 /2p 2 )/(2^ 2 ) 



M(t,Q)^2p 2 /(TTAQ 4 



(8) 



(9) 



U = J sin^ 

fV t 



(g(t, r + R) - g(t, r ) 



d 3 r 
A 3 



, (6) 
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V (x + C/2 - v t , y + Y, z a + Z) d(, (7) 



u t 



N = nX 3 , a = Z e 2 /hv , r = (x ,yo,z ), V a {r ) 
1 r for |r | < A and Vq = for |r | > A. In the re- 
integration, the domain min((j/ + Y) 2 + (zo + Z) 2 , y 2 + 
z o) < r ci = i^Q e2 / m e v o) 2 i s excluded for reasons dis- 
cussed above. The ion function Fi has the same struc- 
ture as F e and is simply obtained by substituting ion 
parameters. M is the convolution of M e and Mi, where 

M eii (t,Q) = J exp(iQ-R/h)F e ^(t,Tl)d 2 R/(2TTh) 2 

are the transverse momentum distributions due to the 
electron-projectile and ion-projectile interaction. In the 
following, most of the discussion is restricted to M e . 

Let us discuss the structure of Eqs. ©,10, which are 
presented here for the first time. The detailed analy- 
sis is quite intricate, and here we give only the main 
results without derivation. We observe that only small 
enough f e can contribute to M and that therefore, ow- 
ing to the large multiplier N in Eq. © , the sin 2 -term 
needs to be small and can be expanded. Then M de- 
pends essentially on the parameter combination a 2 N 
only; more rigorous analysis gives 7 = a 2 A^lniV. The 



for 2pq\uA < Q 2 < (2m e vo) 2 , where A is a solu- 
tion of A = ln(27rainAA 2 wot) ^> 1, ai = min(l,a) 
and Pq = vht. Physically, the Gaussian distribution 
at small Q corresponds to small angle scattering and 
the power-law tail to close collisions with large mo- 
mentum transfer. The tail is obtained only as long as 
t < t\ = 2m 2 e v 2 jvk\YL A, such that 2pglnA < (2m e v ) 2 . 
For longer times, small angle scattering dominates both 
small and large Q regions and the tail disappears. 

We have checked the occurrence of this power-law tail 
by molecular dynamics (MD) simulations. We consider 
a test particle with vq 3> Vth scattered completely classi- 
cally in a finite plasma volume, having dimensions I of the 
order of the screening length. The simulation has been 
performed for a model case with N = 80 and a = 1.9, 
just feasible on a modern PC. Results are plotted in Fig.[21 
for time t — 21/vq < t\ . The histogram presenting the 
MD results is in best agreement with the present the- 
ory (solid curve), clearly showing the power-law tail at 
high momenta. For comparison, also the purely Gaussian 
distribution obtained from the Landau collision integral 
is given as dashed line. Details of these simulations are 
outlined in the caption. Here we should make it clear 
that the power law tail originating from close collisions 
is obtained in nearly identical form within the binary 
collision approach, as it was shown by Landau @ and 
Vavilov 9]. The present theory differs for small-angle 
scattering and therefore in the Gaussian part of M(t, Q). 
To show the difference quantitatively, we have also solved 
the kinetic equation used in 0,0. The result can be writ- 
ten in a form equivalent to Eq. JSJ) with the function f e 

now given approximately by = vt ln(aiA/i?)(R/ft) 2 , 
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C(t) 



;t/(m e v ) 2 J a{Q)d 2 Q and a(Q) 




200 400 600 800 
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FIG. 2: Comparison between MD simulation (histogram), 
present theory (solid line), and predictions of the traditional 
diffusion approximation (dashed line); nM(t,Q) is plotted 
versus Q 2 /ql for time t = 21/ Do, I = 7.239 x 10" 6 cm, velocity 
of test particle vo = 2.297 x 10 8 cm/s, and q\ = l-KnZ^e^t/v®. 
The insert shows the same plot, but with linear scale and 
zoomed to low Q. The simulation assumes an equal num- 
ber of randomly distributed, fixed Coulomb centers of oppo- 
site charge ±e and densities n+ = n_ — n — 1.054 x 10 17 
cm -3 ; the cold plasma limit is chosen with thermal velocity 
vth <S vo. The plasma volume is taken as V — 61 x 21 x 21 
with the test particle (Zo = 2, mo = m e ) moving along 
the central axis in x-direction and starting at a distance 21 
from the surface. The trajectory of the test particle, in- 
teracting with all Coulomb centers, is obtained by solving 
the classical equation of motion by a second-order scheme 
with an adaptive time step. The histogram corresponds to 
2.06 x 10° independent trials. The solid line has been obtained 
numerically from Eqs. JBJ— (7[ with Fi = F e = exp(— / e ). 
The screening length is set to A = I such that the finite 
plasma volume seen by the test particle in this model sim- 
ulation just mimics the physically screened volume occur- 
ing in reality. The straight dashed line is the Landau colli- 
sion integral prediction M(t,Q) = exp(— Q 2 /Qo)/( 7r Qo) with 
Qo = 87re 4 Zo( n + + to—)Lot/vo; Lq — 6.5 is the classical 
Coulomb logarithm evaluated for the parameters of the sim- 
ulation. 



where R = max(o!|R|, ah/m e Vo) and ai = min(l,a). 
The effect of the present theory is that the Gaussian part 
grows more rapidly. This is consistent with enhanced 
small-angle scattering due to simultaneous interaction 
with many plasma particles. 

We have seen in Fig.^ that the quantum limit (7 < 
1) is relevant only in a marginal parameter region. 
Nevertheless, it is contained in M e . For 7 < 1, 
one can use first-order expansions of F e = 1 — f e 
and of the sin 2 -term in f e to find, after some alge- 
bra, M e (t, Q) = C(t)8(Q) + a(Q)nvt/(m e v) 2 , where 



r 2 cl (2m e vo) /{Q 2 + (fi/A) 2 ) 2 is the cross-section of the 
screened Coulomb potential. This leads to (Q 2 ) e — 
AvtL q and the quantum (Bethe) logarithm L q = 
hi(\m e V{)/h). This first-order Born result is obtained 
here for aN 1 / 2 << 1, but not for a < 1 in general. This 
may be understood qualitatively looking at second-order 
processes. Consider the perturbation of the interaction 
of the test particle with a plasma particle j by another 
plasma particle k. This effect is small of order a 2 , but 
since for a plasma with long-range Coulomb forces N 
particles contribute to this second-order process, it can 
be neglected only if a 2 N < 1. 



H, 



Let us finally calculate the energy £(t) = (tp 

the plasma gains due to the energy loss of the test parti 
cle. Here 



Hn 



-h 2 A J /(2m, 




is the Hamiltonian of the plasma without the test par- 
ticle and tp(t, ro, ri, . . . ) the full wavefunction. Mak- 
ing use of the same approximations as in the derivation 
of Eq. JSJl, straightforward algebra leads to d£/dt = 
(Q 2 ) e /2m e + (Q 2 ) i /2mi, where the first term is the 
contribution from plasma electrons with mass m e and 
the second from ions with mass mj. The corresponding 
electron part of the stopping power is then found in the 
standard form 



dE 
dx 



1 d£ 

vq dt 



4irne 4 Z 2 
2 — ' 



but now with L = L c i for 7 > 1 and L = L q for 7 < 1. 

In conclusion, it has been shown that the theory 
of Coulomb scattering in dilute plasma needs revision. 
Bohr's classical Coulomb logarithm L c i is found to apply 
for a^/N > 1 rather than a > 1, and this covers most 
of the density-temperature plane, relevant to practical 
applications. This result calls for experimental verifica- 
tion. We propose to measure energy loss of fully stripped 
ions in carefully characterized, fully ionized plasma lay- 
ers. The parametrically different dependence of L c \ and 
L q on ion charge Zq and velocity vq should allow for a 
clear distinction. 
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